How Much Phase Coherence Does a Pseudogap Need? 
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It has been suggested that the "pseudogap" regime in cuprate superconductors, extending up to 
hudreds of degrees into the normal phase, reflects an incoherent d-wave pairing, with local super- 
conducting order coherent over a finite length scale £, insufficient to establish superconductivity. 
We calculate the single-particle spectral density in such a state from a minimal phenomenological 
disordered BCS model. When the phase-coherence length exceeds the Cooper pair size, a clear 
pseudogap appears. The pseudogap regime, however, is found only over a relatively narrow range 
of phase stiffnesses, hence is not expected to extend more than about 20% above T c . 



It is widely believed that the peculiar normal pseudo- 
gap regime in underdoped cuprate superconductors (see 
rcfs. ||, H for recent reviews), holds keys to unraveling 
the entire problem of cuprate superconductivity. In this 
regime, between the superconducting transition temper- 
ature T c , and T*, possibly hundreds of degrees higher, 
spectral density near the putative Fermi surface is su- 
pressed - in contrast with the familiar behavior of BCS 
superconductors. Evidence of this phenomenon is con- 
sistently observed most clearly in angle-resolved pho- 
toemission (ARPES) , || |[ but also in c-axis tunneling, 
magnetic susceptibility, heat capacity, Raman scattering, 
neutron scattering, and NMR measurements. The vari- 
ation of the pseudogap with momentum, strongest near 
the (ir, 0) directions and weak or nonexistent near (tt, 7r), 
mirrors that of the full d x i_ y i superconducting gap. 

One potential explanation of this behavior, as sug- 
gested by Kivelson and Emery H, is that local supercon- 
ducting correlations (Cooper pairing) set in at T*, but 
that long-range phase coherence is not established un- 
til the temperature drops below T c . The plausibility of 
the idea derives from the facts that phase fluctuations 
generally play a larger role in lower- (here two) dimen- 
sional systems, and the superfluid density in the cuprates 
is very low, making the order parameter phase "floppy." 
The notion is analogous to a magnetic material in which 
local moments form far above the temperature at which 
they become ordered. Nevertheless, this is not quite the 
same thing as pre-formed pairs. 

Several previous calculations have investigated this 
idea[|[ 0, H[ |). With the exception of reference fj], which 
studies effects of uniform superflow fluctuations, these 
are essentially diagrammatic approaches, with the atten- 
dent need for various sorts of truncations or resumma- 
tions. The purpose of the present work is to get a better 
understanding of the possibilities and limitations of the 
incoherent pairing scenario in a BCS-like framework with 
as few additional uncontrolled approximations or dubious 
assumptions as possible. We do that by investigating a 
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FIG. 1: spectral weight at ((3/24)tt, tt) for A = 0.05t. In- 
dividual energy levels are broadened by convolution with 
a Lorentzian of FWHM t/100. The correlation length is 
computed from the Kosterlitz-Thouless expression £/a = 
0.15 exp (l.82\/t) @, 0, with t = K/(l - K). 



minimal phenomenological model permitting essentially 
exact (albeit numerical and statistical) solution. This 
model has the following ingredients. 
(1) Electrons propagate according to a tight-binding 
model in the presence of a disordered pair order param- 
eter with local d x 2_ y 2 symmetry: 
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We take the ratio of 



next-nearest-neighbor hopping (t') to nearest-neighbor 
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(t) to be t'/t = — 0.35|l(J, and the chemical potential 
/i = — 1.1076t, corresponding to a doping of x = 0.15. 
We estimate t r* 1/5 eV, roughly 1/4 of the observed 
bandwidth. 

(2) The complex pair field Ay (= Ajj) couples to singlet 
pairs spanning nearest-neighbor sites i and j. In a d x 2_ y 2 
superconductor, the (spatially uniform) pair field Ay is 
Ao on horizontal bonds, and — Ao on vertical bonds. 

(3) In a traditional BCS approach, the gap parameter 
Ao is computed self-consistently taking into account the 
underlying pairing interaction. In our phenomenological 
approach, this linkage is decoupled. Instead, the pair 
field Ay = A e l6,ij in eq. (|l|) is treated as a quenched 
random variable, chosen from a thermal ensemble. We 
assume that |Ay |, the magnitude of the pairing correla- 
tions is a fixed (large) energy scale Ao equal to the full 
gap below T c , but that the phase suffers thermal fluctu- 
ations determined by the classical XY model 

fiHxY = 1.12 K 

cos(0ij - Oik). (2) 

(ij,ik) 

The sum runs over all pairs of nearest-neighbor links, i.e., 
links ij and ik related to each other by a 90° rotation 
about a common site i. The centers of these links form 
a square lattice, canted at 45° to the lattice of tight- 
binding sites. The factor of 1.12 in eq. (||) means that 
K is measured in units of the critical coupling, i.e., the 
Kosterlitz-Thouless-Berezinskii (KTB) transition occurs 
at K = 1 |nj . One may think of the Hamiltonian of Eq. 
^| as arising from the Hamiltonian 

Ti-pair = 9 ^ Ay- Ajfc, (3) 
i;jk 

which describes the pivoting |l3| of a Cooper pair from a 
nearest-neighbor link ik to a perpendicular link ij (see 
inset of Fig. ||). Positive g encourages local d- wave order. 

Assumptions 1-3 above specify our model. In con- 
trast to BCS theory, a distinction is drawn between 
the local pairing strength Ao and the phase stiffness 
K; we take them as independent phenomenological pa- 
rameters. A genuinely two-dimensional system of this 
sort develops only algebraic order for K > 1. In a 
real material, coupling between planes causes crossover 
to three-dimensional behavior and true long-range or- 
der. For K < 1, there is no superconductivity, but local 
superconductor-like correlations. The primary questions 
are: for what range of phase stiffness K does the local 
pairing manifests itself above the KTB transition, and, 
to how high a temperature is the picture able to hold 
together? 

Methods. Given a pair field Ay , we solve numerically 
for the ground state and quasiparticles by a suitable Bo- 
goliubov transformation: 
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FIG. 2: Upper right-hand quadrant of the Brillouin zone grid 
used for numerical computation. The heavy solid curve in- 
dicates the normal state Fermi surface for the parameters 
(t,t',/j,) reported in the text. The dashed curves are normal- 
state energy contour lines at a spacing of t/2. Circled points 
are those for which spectral density results are displayed in 
this paper. Inset: a pair hopping event. 

where the index i labelling lattice sites is summed over on 
the right-hand side. Since the T operators destroy quasi- 
particles, we have the Bogoliubov-de Gennes equation 

where and are column vectors to be determined. 

With Bogoliubov quasiparticle states in hand, we 
calculate the experimentally relevant particle and hole 
spectral densities A + (kw) = (c(k)8(H — u>)c(ky) and 
A~(kui) = (c(k)5(H — uj)c{k)^). Since the basic as- 
sumption is that the crucial thermal effects are fluctu- 
ations of the order parameter phase, we compute the 
quenched average spectral densities at zero "quasiparticle 
temperature" H| over the distribution of {Ay}: 




Representative full spectral densities Ao(k,oj) [equal 
to A^(k,— lo) for uj < and to A^{k,uj) for u > 0] 
over a range of K are displayed in Figures [l] and ^ for 
a point near the Fermi surface in the (n, 0) direction 
with Ao/t = 0.05 and 0.15, respectively. Results for a 
Fermi surface point much closer to the (jr, w) (nodal) di- 
rection is shown in Figure ^| These plots are the result 
of Monte Carlo sampling of the pairing field ensemble 
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FIG. 3: spectral weight at ((3/24)7T, it) for A = 0.15i. 

for a 48 x 48 site lattice, with averaging over three to 
five disorder realizations and crystallographically equiv- 
alent k points. In each figure, the inner dashed lines 
indicate plus or minus the tight-binding energy e(k), 
and the outer dashed lines are the quasiparticle ener- 
gies E (k) = ±y / e(k) 2 + A (fc) 2 for a perfectly ordered 
pairing field with the specified Ao (i.e., K — oo). 

What do these results tell us? First, the single-particle 
spectral function exhibits a pseudogap even for fairly 
short-ranged phase coherence. Roughly speaking, the 
phase correlation length, £, is the distance over which 
Cooper pairs propagate coherently, or the characteristic 
vortex-antivortex separation. The figures show a pseu- 
dogap developing by about £ = 15a for Ao = 0.05t. For 
this pairing magnitude, the Cooper pair size, estimated 
according to £o = fcv F /A (k), is of order 8a. (For this 
estimate, we use Fermi velocity vp and full gap at the 
Fermi surface near the (ir, 0) direction. Other choices 
give an estimate a few times larger.) 

Second, the onset of the pseudogap with K is rela- 
tively rapid near the KTB critical point. The (ir, 0) gap 
achieves nearly half its full value by the KTB transi- 
tion, continuing to increase slowly with increasing K. 
The variation of the pseudogap along the Fermi sur- 
face closely resembles the pure d x 2_ y 2 form. If we in- 
vert the usual BCS relationship between energy gap 
and pairing amplitude, using the peak position of the 
spectral density as a surrogate for the BCS quasiparti- 
cle energy, we can extract an effective "gap parameter" 
A eff (fc) = v/Spca^fc) 2 -e(fc) 2 . At K = 0.9, and for 
both A /t = 0.05 and A /t = 0.10, the ratio A off /A (fc) 
of effective gap to full gap is between 0.41 and 0.63 ev- 
erywhere on the Fermi surface (with most points being 



FIG. 4: Spectral weight at (7, 13) fj for Ao = 0.05i. Note the 
curious feature that the peak is below the A = value for 
smaller values of K. The average, however, is closer because 
of the long high-energy tail. 



close to 0.5). 

Third, spectral peaks are broad with widths roughly 
linear in the pairing magnitude Ao- As shown in Figure 
|3|, this broadening can wash out the clear pseudogap, 
though it does not affect the shift of the peaks. The 
widths do not monotonically decrease as K is increased, 
but are greatest near K = 1. 

Fourth, although the pseudogap appears at a modest 
phase correlation length (comparable to a few Cooper 
pair radii), the most natural interpretation of this in 
terms of temperature is quite discouraging if the aim 
is to explain experimentally observed pseudogaps. As 
stressed earlier, K is a phenomcnological parameter, but 
it should probably be close to \A \ 2 n®/m, where n° is 
the bare superfluid density. In that case, K is roughly 
proportional to 1/T, and our results imply T*/T c < 1.2. 
The collapse of the pseudogap in our model corresponds 
to £ ~ £o = ^w F /Ao, the Cooper pair size. If a psuedogap 
is to be maintained to high temperature, something must 
halt the decrease of £ at that point. In other words, the 
order parameter phase must remain stiff at wavelengths 
comparable to £o- This issue is not addressed by our 
phenomcnological model. 

It is interesting to see how well the numerical re- 
sults can be reproduced by a standard diagrammatic 
approximation, using a self-consistent Green function. 
This approximation is similar to that used by Kwon and 
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DorseyH, and is represented by 



(^=)" 1 =M" 1 + 



,-X-.. 



Equivalently, the self-energy is given by 

nG](k,u) = [(\A(q)\ 2 )G(-u;,q-k)d 2 q. (7) 

We solve this iteratively, including (|A(0)| 2 ) only in 
the final iteration. For K < 1, the pair field cor- 
relator is accurately approximated by (|A(q)| 2 ) oc 
[£~ 2 + 2 X) e (l — cos(q • e))] , with a proportionality con- 
stant fixed by the sum rule J2 | A(q)| 2 = 27V A 2 , (N is the 
number of sites in the lattice). Spectral densities calcu- 
lated via this method (Fig. ||) also show the development 
of a pseudogap at values of K in close agreement with 
the exact results. However, the peaks are narrower, by a 
factor of as much as two. 

In summary, we have shown, using a simple phe- 
nomenological model, that a single-particle pseudogap, 
per se, is natural in the presence of order parameter 
phase fluctuations. In this context, it is unclear what 
mechanism could maintain short range phase stiffness 
up to the experimentally observed T*. On the other 
hand, recent microwave conductivity measurements! 15 
suggest that the bare superfluid density vanishes at a 
temperature comparable to the most natural pseudogap 
collapse temperature found here, about 20% above T c . 
This is also consistent with other previous theoretical 
findings @| | . 
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FIG. 5: Spectral density at (3,23)|j for A = 0.05t and 
several values of the phase correlation length £ (lattice- 
spacing units), calculated by the self-consistent Green func- 
tion method. 
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